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We further study the "complementary" Ansatz, Tr(M„)=0, for a prediagonal light Majorana type 
neutrino mass matrix. Previously, this was studied for the CP conserving case and the case where 
the two Majorana type CP violating phases were present but the Dirac type CP violating phase was 
' neglected. Here we employ a simple geometric algorithm which enables us to "solve" the Ansatz 

, including all three CP violating phases. Specifically, given the known neutrino oscillation data and 

/"S ■ an assumed two parameter (the third neutrino mass ms and the Dirac CP phase 5) family of inputs 

we predict the neutrino masses and Majorana CP phases. Despite the two parameter ambiguity, 
interesting statements emerge. There is a characteristic pattern of interconnected masses and CP 
phases. For large 7713 the three neutrinos are approximately degenerate. The only possibility for a 
mass hierarchy is to have 1TI3 smaller than the other two. A hierarchy with 7713 largest is not allowed. 
Small CP violation is possible only near two special values of 777,3. Also, the neutrinoless double 
beta decay parameter is approximately bounded as 0.020 eV < |m ee | < 0.185 eV. As a byproduct 
of looking at physical amplitudes we discuss an alternative parameterization of the lepton mixing 
matrix which results in simpler formulas. The physical meaning of this parameterization is explained. 
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O • I. INTRODUCTION 

^ \ The remarkable experimental achievements (for some recent examples see refs. 0, 0> 0) relating to neutrino 
. oscillations QJ have brought much closer to reality the goal of determining the "light" neutrino masses and the 
presumed 3x3 lepton mixing matrix. It is possible that more than three light neutrinos are required in order to 
understand the results of the LSND experiment j5j. However, we consider it reasonable, before deciding on this, to wait 
for further supporting evidence as should be supplied soon by the MiniBooNE collaboration The mixing matrix 
contains three mixing angles and, if the neutrinos are considered to be Dirac type fermions, a single CP violation 
phase. That would be completely analogous to the situation prevailing in the quark sector of the electroweak theory. 
But it seems very interesting to consider the possibility that the three light neutrinos are Majorana type fermions. 
This involves only half as many fermionic degrees of freedom and would be mandated if neutrinoless double beta 
decay were to be conclusively established. The Majorana neutrino scenario implies the existence of two additional 
CP violation phases[3,ll,H,ll3- Then nine quantities (beyond the charged lepton masses) would be required for the 
specification of the lepton sector: three neutrino masses, three mixing angles and three CP violation phases. 

According to a recent analysis [Tl| it is possible to extract from the data to good accuracy, two squared neutrino 
mass differences: m| — m\ and |rn§ — m\\, and two inter-generational mixing angle squared sines: s\ 2 and s^- 
Furthermore the inter-generational mixing parameter s\ s is found to be very small. Thus 5 out of 9 quantities needed 
to describe the leptonic sector in the Majorana neutrino scenario can be considered as "known" . For many purposes 
it is desirable to get an idea of the remaining 4 parameters. As an aid in partially determining the other parameters, 
a so-called "complementary Ansatz" was proposed E^> Q El ■ The name arises from the fact that if CP violation 
is neglected, the Ansatz determines (up to two different cases) all three neutrino masses, given the two known squared 
mass differences. 
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This complementary Ansatz simply reads, 

Tr(M v ) = 0. (1) 

Here M„ is the symmetric, but in general complex, prediagonal Majorana neutrino mass matrix. It is brought to 
diagonal form by the transformation, 

U T M V U = diag(mi,m2,m 3 ) — M v , (2) 

where U is a unitary matrix and the mj may be chosen as real, positive. We impose the condition in a basis where 
the charged leptons are diagonal so that U gets identified with the lepton mixing matrix. 

Since Eq. (JJJ comprises two real equations, it gives two conditions on the unknown 4 parameters of the lepton sector. 
In other words, the lepton sector is described by a two parameter family of solutions. This can be approximately 
simplified 0,0, El by noting that the effects of one CP violation phase get suppressed when the small quantity S13 
vanishes. Then there is a one parameter family of solutions describing the lepton sector and it is straightforward to 
compute physical quantities for parameter choices which span this family. The main purpose of the present paper is 
to find the general solutions of Eq. Q without making this approximation. This gives a two parameter family which 
allows one to study the interplay of all three CP violation phases. 

A plausibility argument supporting the complementary Ansatz is concisely presented in sections 2 and 3 of |Tfjj . It 
is based on using the SO(10) grand unification group in the approximation that the non-seesaw neutrino mass term 
dominates. The Higgs fields which can contribute to fermion masses at tree level are the 10, 120 and the 126. If 
only a single 126 appears (but any number of the others) one has the relation 

Tr(M D - rM E ) oc Tr(M v ), (3) 

where M D and M E are, respectively, the prediagonal mass matrices of the charge -1/3 quarks and charge -1 leptons, 
while r sa 3 takes account of running masses from the grand unified scale to about 1 GeV. Now one of the major 
surprises generated by the neutrino oscillation experiments is that, unlike the quark mixing matrix which has the 
form diag (1,1,1) + 0(e), the lepton mixing matrix is not at all close to the unit matrix. This suggests a further 
approximation in which one takes M D and M E to be diagonal but allows the neutrino mass matrix to be far from 
the unit matrix. Then the left hand side of Eq. @ is approximately (nib — 3m T ), which is in turn close to zero. 

As we will see, the model makes a number of characteristic predictions for the neutrino mass spectrum which should 
enable it to be readily tested in the near future. A very recent review of many other models is given in ref. |17| . 

For convenience, our notation (essentially the standard one) for the lepton mixing matrix and the corresponding 
parameterized Ansatz is given in section II. 

In section III the Ansatz is solved in the sense of providing a geometrical algorithm which, given the two input 
quantities 7773 (third neutrino mass, taken positive) and 5 (conventional CP violation phase in the lepton mixing 
matrix), enables one to predict the other two neutrino masses as well as the other two CP violation phases. Of course, 
the experimental knowledge on the neutrino squared mass differences and CP conserving intergenerational mixing 
angles are taken to be "known" . We separate the solutions into two types I and II, depending respectively on whether 
7773 is the largest or the smallest of the neutrino masses. In addition, there is a discrete ambiguity corresponding to 
reflecting a triangle involved in the algorithm. A "panoramic" view of the predictions as functions of 7773 and 6 are 
presented in a convenient tabular form. The greatest allowed value of 7713 is determined by a cosmology bound. As 
?7i3 decreases, a point is reached at which the type I solutions no longer exist. As 7773 decreases even further, the type 
II solutions also cease to exist. The corresponding values of 7773 at which these solutions become disallowed depend 
on the assumed value of the input 5. This correlation is studied analytically. 

Some physical considerations are discussed in section IV. First, the dependence on the experimentally bounded 
squared mixing angle, s^ 3 is investigated. We present also a chart showing the dependence of the neutrinoless double 
beta decay parameter \m ee \ on the input parameters 7773 and 5. Even though the inputs are varying over a fairly large 
range, the rather restrictive approximate bounded range 0.020 eV < |m ee | < 0.185 eV emerges from the Ansatz. 

After calculating observable quantities in the model one observes that they depend more simply on certain linear 
combinations of the conventional "Dirac" and "Majorana" CP violation phases. In section V we discuss an alternative 
parameterization of the lepton mixing matrix in which these combinations occur directly. In this parameterization 
the three phases just correspond to the three possible intergenerational mixings. An "invariant" combination of these 
three corresponds to the usual "Dirac" phase 5. 

We conclude in section VI which contains a brief summary and a discussion of results which emphasize some unique 
features of the present work. 
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II. PARAMETERIZED COMPLEMENTARY ANSATZ 

We define the lepton mixing matrix, K from the charged gauge boson interaction term in the leptonic sector of the 
electroweak Lagrangian: 



C=^W-e Llli Kv + H.c. 



(4) 



Note that the "mass diagonal" neutrino fields, are related to the fields, pi in the prediagonal mass basis by the 
matrix equation, 



p = Uv. 

We adopt essentially what seems to be the most common parameterization: 

K = KexpLUQ 1 (t), 

where a unimodular diagonal matrix of phases is defined as, 

«o(r) = diag(e iT \e^,e tT3 ), 

n+ r 2 +T3 = 0. 



(5) 
(6) 

(7) 



The remaining factor, K exp which is the only part needed for describing ordinary neutrino oscillations is written as 
the product of three successive two dimensional unitary transformations, 



K exp = ^23(^23, O)wi 3 (0i 3 , -£)wi 2 (0i2, 0), 
with three mixing angles and the CP violation phase 6. For example in the (12) subspace one has: 



wi 2 (0i2, ^12) 



cos6 12 e^ 12 sin0 12 



(8) 



(9) 



1 

with clear generalization to the (23) and (13) transformations. Multiplying out yields: 



Kr 



where = sinOij and Cy = cost 



C12C13 



S12C13 



«i3e 



-i<5 



-S12C23 - ci 2 Si 3 s 2 3e lS C12C23 - Si 2 si 3 s 23 e lS ci 3 s 23 



S12S23 - Cl2Sl3C23e 



IS 



-C12S23 - Sl2Sl3C23e l<5 C13C23 



(10) 



Identifying K with U in Eq. the Ansatz of Eq. now reads, 

Tr(M v K-} p K; xp LO {2 T )) = 0, 



(11) 



where Eqs. J2J and © were used. With the parameterized mixing matrix of Eq. (|10l) the Ansatz finally becomes: 



rti\e 



2iri 



[1 - 2i(c 12 s 13 ) 2 sin6e- lS ] + 
m 2 e 2lT2 [l - 2i(s 12 s 13 ) 2 sinSe" lS ] + 
m 3 e 2tT3 [l + 2i{s 13 ) 2 sinSe tS ] = 0. 



(12) 



In this equation we can choose the diagonal masses mi,m 2 ,m 3 to be real positive. Notice that setting the mixing 
parameter S13 to zero eliminates the dependence on the CP violation phase 5. Then Eq. I|12fl goes over to the simpler 
form studied previously 0] . 



III. SOLVING THE ANSATZ IN THE GENERAL CASE 

For dcfinitcncss we will use the following best fit values for the differences of squared neutrino masses obtained in 
ref. [13: 

A = m\-m\ = 6.9 x 10" 5 el/ 2 , 
B=\m\-m 2 2 \ = 2.6 x 10^ 3 ey 2 . (13) 



4 



The uncertainty in these determinations is roughly 25%. Similarly for definiteness we will adopt the best fit values 
for a\ 2 and s 23 obtained in the same analysis: 

4, = 0.30, 4=0.50. (14) 

These mixing angles also have about a 25% uncertainty. The experimental status of s 2 3 is less accurately known. At 
present only the 3 a bound, 

S ? 3 < 0.047, (15) 

is available. For our discussion we will consider s 2 3 to be known at a "typical" value satisfying this bound and examine 
the sensitivit y to changing it. Of course, the experimental determination of s±3 is a topic of great current interest. 

Previously |l5j, the (positive) mass of the third neutrino, 7713 was considered to be the free parameter. It was varied 
to obtain, via the simplified Ansatz equation, a "panoramic view" of the two independent Majorana phases (say T\ 
and t 2 ). In the present case, we shall not neglect the CP violation phase 5 and consider it too as a free parameter 
to be varied. It is necessary to specify a suitable algorithm to treat the full Ansatz. Previously it was noted that 
the simplified Ansatz could be pictured as a vector triangle in the complex plane having sides equal to corresponding 
neutrino masses (See Fig. 1 of [l5|). The three internal angles were found by trigonometry and related to the three 
angles made by the sides with respect to the positive real axis. Those in turn were twice the three (constrained) 
Majorana phases. The orientation of the triangle got determined (up to a reflection) by the constraint in Eq. (J7J. In 
the present case we will also rewrite the Ansatz equation as a vector triangle in the complex plane. However the sides 
will differ from the neutrino masses. In addition the angles will differ from twice the constrained Majorana phases. 

To start, we choose a value for 7713 and a value for the phase 5. Then we can obtain from Eas. H13(l two different 
solutions for the other masses toi and m 2 . We call the solution where TO3 is the largest neutrino mass, the type I 
case. The case where TO3 is the smallest neutrino mass is designated type II. m\ will be determined from the assumed 
value of 7713 as, 

m\ = m\ - A =F B, (16) 
where the upper and lower sign choices respectively refer to the type I and type II cases. In either case we find m 2 as, 

m\ = A + m\. (17) 

Next, we redefine variables so that each of the three terms in the Ansatz equation, i|12|) is characterized by a single 
magnitude, to- and a single phase, 2t[. Eq. (|T2t then reads, 

m\e 2i ^ + m' 2 e 2iT * + m^e 2 "^ = 0. (18) 
This equation evidently represents a vector triangle in the complex plane, as illustrated in Fig. ^ However the lengths 




FIG. 1: Vector triangle representing Eq. 118t . 



are not the physical neutrino masses and the phases are not twice the physical Majorana CP violation phases. The 
auxiliary, primed, masses are seen to be related to the (now known) physical masses by, 



TO- = GiTOi, 



(19) 
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l/ 2 

1-4 (ci 2 si 3 ) z sin"* 5 + 4 (c^Sis)^ sin 2 (5 . 



where, 

Gi = 

G 2 - 

G 3 = [l-4s2 3 sin 2 5 + 44 3 sin 2 (5] 1/z . (20) 



4 (S12S13) 2 sin 2 5 + 4 (s 12 si 3 ) sin 2 5 



1/2 



Notice that, since 6 has already been specified, the relations between the m! i and the rrii are now known. Similarly, 
the physical phases are related to the primed ones by, 



Fi, (21) 



where, 



1 (ci2Si 3 ) 2 sin(2(5) 



F\ = — arctan 2 , 

2 1-2 (C12S13) sin <5 

1 (s 12 si 3 ) 2 sin(2<T) 

r2 — — arctan ■ 



2 ' l-2(s 12 si 3 ) sm 2 <5 

1 s? 3 sin(2(5) 
F 3 = — arctan — — 5-. (22) 
3 2 l + 2s 2 3 sin 2 <5 V ' 

Again, note that, since 5 has been specified, the relations between the and the r 8 ' are now known. Referring to Fig. 
1 we can determine the internal angles, p[ by using the law of cosines. For example, 

cos ^ = 2^4 ' (23) 

Next, the auxiliary phases t[ can be related to the internal angles just obtained as: 

t[ = ^(tt-Mi - 2 ^ 2 ) +p, 
t! 2 = i(7r + 2// 1 +/4)+p, 

4 = i(-27r-/ii+/4)+P- (24) 

The remaining still unknown parameter here is p which we added to the right hand side of each equation. It represents 
the effect of an arbitrary rotation of the whole triangle, which should not be determinable from the internal angles. 
It can be determined, however, by making use of the constraint on the physical phases, Tj = 0. Notice that there 
is no corresponding constraint for t[. Using Eq. (|21|l we get, 

i 

where the Fi are to be read from Eqs. I|22|l. Now the masses mi, m2 and the phases ti,T2 have been determined by 
a simple algorithm, upon specification of m 3 and <5. 

As remarked above, the dependence on the CP violation phase 8 is suppressed in the limit that the mixing parameter 
sf 3 vanishes. Hence, to illustrate this new feature, we will consider a value s 2 3 = 0.04, close to the 3<r upper bound of 
0.047 llj. The predictions of the neutrino masses (mi, m 2 ) and two independent phases (tj, r 2 ), from the Ansatz for 
various assumed values of m 3 and 6 are given in Table [I] Representative values of S were chosen to lie between and 
7r since it may be observed from Eqs. I|20() and 1L'2I) that the solutions will have a periodicity of tt with respect to S. 
Just from the Ansatz there is no upper bound on the value of m 3 . However there is a recent cosmological bound 0] 
which requires, 

\mi\ + \m 2 \ + |ms| < 0.7 eV. (26) 

Thus values of m 3 greater than about 0.3 eV are physically disfavored. Table |U shows that at this value both type I 
and type II solutions exist. This is true also for higher values of rn 3 . The picture remains very similar down to around 
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type mi, 7712,7713 in eV n, T2(S = 0) ri, T2(<5 = 0.5) n, 7-2(6 = 1.0) Ti,T2(5 = 1.5) Ti,T2(5 = 2.0) n, T2(c5 = 2.5) 



I 


0.2955, 


0.2956, 


II 


0.3042, 


0.3043, 


I 


0.0856, 


0.0860, 


II 


0.1119, 


0.1123, 


I 


0.0305, 


0.0316, 


II 


0.0783, 


0.0787, 


II 


0.0643, 


0.0648, 


II 


0.0541, 


0.0548, 


II 


0.0506, 


0.0512, 


II 


0.0503, 


0.0510, 



0.3 

0.3 

0.1 

0.1 

0.06 

0.06 

0.04 

0.02 

0.005 

0.001 



0.0043 

-0.0041 

0.0486, 

-0.0311 

0.3913, 

-0.0669 

-0.1064 

-0.1747 

-0.2601 

-0.3830 



1.0428 
1.0512 
0.9975 
1.0774 
0.6543 
1.1119 
1.1494 
1.2115 
1.2620 
1.1805 



0.0126 

0.0043 

0.0566, 

-0.0226 

0.3873, 

-0.0583 

-0.0978 

-0.1669 

-0.2603 



1.0495 
1.0577 
1.0049 
1.0835 
0.6748 
1.1174 
1.1541 
1.2142 
1.2611 



0.0058, 
-0.0023 
0.0489, 
-0.0289 
0.3578, 
-0.0644 
-0.1040 
-0.1751 
-0.2914 



1.0536 
1.0615 

1.0106 
1.0863 

0.7048 

1.1188 
1.1538 
1.2095 
1.2276 



-0.0108 
-0.0189 
0.0318, 
-0.0453 
0.3288, 
-0.0806 
-0.1203 
-0.1928 
-0.3251 



1.0510 
1.0587 

1.0091 
1.0828 

0.7167 
1.1145 
1.1483 
1.2012 
1.2035 



-0.0210, 1.0440 

-0.0291,1.0518 

0.0219,1.0015 

-0.0556,1.0763 

0.3258,0.7013 

-0.0911,1.1085 

-0.1307,1.1430 

-0.2024,1.1976 

-0.3250,1.2094 



-0.0153, 1.0394 

-0.0235,1.0476 

0.0285,0.9953 

-0.0503,1.0731 

0.3530,0.6720 

-0.0860,1.1066 

-0.1255,1.1428 

-0.1951,1.2019 

-0.2950,1.2369 



TABLE I: Panorama of solutions as 7713 is lowered from about the highest value which is experimentally reasonable to about 
the lowest value imposed by the model. For each value of 7773, predictions are given for (7711,7712) and for (n, T2) in the cases 
where <5 = 0,0.5, 1.0, 1.5,2.0,2.5. All phases are measured in radians. Here the value sf 3 =0.04 was adopted. In the type I 
solutions 7773 is the largest mass while in the type II solutions 7773 is the smallest mass. 



777.3 = 0.1 eV but as one gets closer to roughly 0.06 eV, there is a marked change. If one further lowers 7773, it is 
found that the type I solution no longer exists. On the other hand the type II solution persists and does not change 
much until 7773 approaches the neighborhood of 0.001 eV. There are no solutions for 777,3 below this region. 

Note that the columns in Table |U with 6 = correspond to the previous case, discussed in some detail in section 
IV of 15]. In this case, m! i and r[ respectively coincide with m, and so we can identify the vectors of the triangle 
with the physical masses and phases. As one decreases the value of 7773 the type I triangle goes from being close 
to equilateral to the degenerate situation with three collinear vectors. In this limiting case the vectors representing 
neutrino one and neutrino two are approximately equal and add up to exactly cancel the vector representing neutrino 
3. The precise orientation of the straight line is due to imposing the constraint in Eq. J7J). This is actually a CP 
conserving case Then one can find the 7773 value (a little below 0.06 eV) for this situation by looking for a real 
solution of mi + 777 2 + 7773 = together with Eqs. ifHfl) (See Eq. (4.4) of pjj). Clearly there can be no type I solutions 
below this value of 7773. The type II solutions can exist below this value but similarly end (a little below 7713 =0.001 
eV) when the triangle becomes degenerate in a different way. For the type II degenerate triangle, the neutrino 1 and 
neutrino 2 vectors are collinear but oppositely directed and the small neutrino 3 vector adds to the neutrino 1 vector 
to cancel the neutrino 2 vector. This is also a CP conserving case. 

When the effects of 5 not equal to zero are included, it is not possible to make a triangle out of the physical neutrino 
masses and phases. The relevant auxiliary triangle is made, as illustrated, using the primed masses and phases. Thus 
the limiting values of 7773, where the type I and type II cases each end, correspond to this primed triangle becoming 
degenerate. We can get the limiting value by looking for real solutions of J2i Gjmj=0, together with Eqs. (|T3|) . The 
limiting value, {m^) m i n is found to be: 

(™ 3 )L n -7^[-/3-(/3 2 -4a 7 ) 1/2 ], (27) 

where, 

a = {Gl + Gl-G%f-AG\Gl 

13 = -2(G 2 + G 2 -G 2 ){AG 2 ±B{G 2 + G 2 2 ))+4G 2 G 2 (A±2B) 7 

7 = (AG 2 ±B(G 2 + G 2 )) 2 t4G 2 G 2 B(A±B). (28) 

Here the upper and lower sign choices respectively refer to the type I and type II cases. 

The computed values of (rri3) m i n as a function of <5 are shown in Table [n] Looking at this table, one can see why 
the entries in Table[I]for 7773 =0.001 eV and non-zero values of 5 are missing. Simply, for those cases, (777,3)™™ > 0.001 
eV. Clearly, this correlation of the allowed input values of 777,3 with the input values of S must be respected in studying 
the present model. This correlation is imposed by the Ansatz itself. For large values of 7713 there is no constraint 
from the Ansatz but Eq. (|26|l gives an experimental constraint. It should be remarked that the collinear auxiliary 
triangles with non-zero S correspond to CP violation. In these cases the CP violation arises from non trivial phases, 
Ti in addition to the assumed non-zero S. Actually, a better measure of CP violation involves the (two independent) 
phase differences r, — Tj. These are the objects essentially related to the internal angles in Fig. ^ 

As noted in ref the possibility of reflecting the triangle about any line in the plane gives another set of 
solutions corresponding to reversing the signs of all the phase differences r, — Tj. In the present case, where S is not 
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type 


{m-i) min (S = 0) in 


eV (m 3 ) m i„ (5 = 


0.5) (m a ) min (8 = 


1.0) (m 3 ) min (5 = 


1.5) (m 3 ) mi „ (5 


= 2.0) {m 3 ) mi n (<5 = 2.5) 


I 


0.0592716 


0.0590967 


0.0587178 


0.0584799 


0.0586203 


0.0589971 


II 


0.0006811 


0.00105461 


0.0019024 


0.0024636 


0.0021294 


0.0012723 



TABLE II: Minimum allowed value of the input mass, 7713 as a function of the input CP violation phase, 5. These correspond 
to the cases where the triangle in Fig. 1 becomes degenerate. Here, the choice S13 = 0.04 has been made. 



zero, reflecting the "unphysical" triangle about any line in the plane will give an alternate solution in which the r[ — rj 
are reversed in sign. More specifically, one should reverse the signs of the first terms on the right hand sides in Eqs. 
(|24|l . The physical phases, Tj for this alternate solution will then depend on 6 as illustrated in Table ITTT1 

type m 3 in gV 71,7-2 (3 = 0) t u t 2 {5 = 0.5) n,-^ = 1.0) r t , 7 2 (<5 = 1.5) n, t 2 (<? = 2.0) n, r 2 (5 = 2.5) ~ 

I 03 -0.0043,-1.0428 0.0113,-1.0393 0.0210, -1.0422 0.0151, -1.0492 -0.0015, -1.0536 -0.0123, -1.0512 

II 0.3 0.0041,-1.0512 0.0196, -1.0475 0.0291, -1.0500 0.0232, -1.0569 0.0066, -1.0614 -0.0040, -1.0593 

I 0.1 -0.0486, -0.9975 -0.0326, -0.9947 -0.0221, -0.9992 -0.0275, -1.0073 -0.0444, -1.0111 -0.0560, -1.0070 

II 0.1 0.0311, -1.0774 0.0465, -1.0732 0.0557, -1.0748 0.0495, -1.0810 0.0331, -1.0859 0.0228, -1.0848 

I 0.06 -0.3913, -0.6543 -0.3634, -0.6645 -0.3310, -0.6934 -0.3246, -0.7149 -0.3483, -0.7109 -0.3806,-0.6837 

II 0.06 0.0669, -1.1119 0.0822, -1.1071 0.0912, -1.1073 0.0849, -1.1127 0.0685, -1.1180 0.0584, -1.1183 
II 0.04 0.1064, -1.1494 0.1217, -1.1438 0.1308, -1.1423 0.1246, -1.1465 0.1082, -1.1526 0.0979, -1.1506 
II 0.02 0.1747, -1.2115 0.1909, -1.2040 0.2019, -1.1981 0.1971, -1.1994 0.1799, -1.2072 0.1676, -1.2136 
II 0.005 0.2601, -1.2620 0.2853, -1.2447 0.3182, -1.2162 0.3294, -1.2017 0.3024, -1.2190 0.2694, -1.2486 
II 0.001 0.3830, -1.1805 



TABLE III: Panorama of solutions, using the reflected triangle, as mz is lowered from about the highest value which is 
experimentally reasonable to about the lowest value imposed by the model. Notice that the predicted masses, mi and m 2 have 
not been given since they are the same as in Table Q] 

Unlike the S = case, the phase differences for the reflected triangle solution, are now only approximately the 
negatives of those for the original solution. For example, in the case of a type I triangle with 7773 =0.3 and 8 =1.0, 
Table|Jshows n — 7"2=-1.0478 for the original solution while Table ITTT1 shows t\ — T2=+1.0632 for the reflected triangle 
solution. 

It should be remarked that the number of decimal places to which we are calculating is chosen in order to be able 
to compare various solutions of the Ansatz with each other for precisely fixed values of the input mass differences and 
mixing angles. The experimental accuracy of the inputs must, of course, be kept in mind. 



IV. PHYSICAL APPLICATIONS 



It is very interesting to note the dependence of our results on the value of the necessarily small quantity s 2 3 , which 
can be seen from the Ansatz Eq. (|12[) to modulate the S dependence. For this purpose let us consider, instead of the 
value 0.04, the value 0.01. The resulting analog of Table^is presented in Table llVl 

Notice that Table I1VI has fewer missing solutions for the case 7773 =0.001 than does Table [I] This is because 
decreasing s 2 3 brings the Gi in Eq. (|19|l closer to unity, which in turn brings the physical neutrino masses closer to 
the auxiliary m\. The modified lower limits for 7713 are illustrated in Table Ivl 

The implications of this model are relevant to experiments which are designed to search for evidence of neutrinoless 
double beta decay. The amplitudes for these processes contain a factor, m ee , which is independent of the nuclear 
wave functions. Its magnitude is given by, 

|m ee | = \ mi (K expll ) 2 e- 2 ^ +m 2 {K expl2 fe- 2 ^ + m 3 (K expl3 ) 2 e - 2 ^\, (29) 

which appears to require, for its evaluation, a full knowledge of the neutrino masses, mixing angles and CP violation 
phases. The present experimental bound 20] on this quantity is 

\m ee \ < (0.35 1.30) eV, (30) 

A very recent review of neutrinoless double beta decay is given in ref . pl| . Using the general parameterization of Eq. 
© one finds. 



m ee = V C 2 + D*, 



(31) 
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type m 3 in eV n,r 2 (8 = 0) ti,t 2 (<S = 0.5) n, t 2 (5 = 1.0) ti,t 2 (5 = 1.5) n, t 2 (<5 = 2.0) n, t 2 (<5 = 2.5) 

I 03 -0.0043,-1.0428 0.0063, 1.0445 0.0046, 1.0455 0.0007, 1.0448 -0.0018, 1.0432 -0.0006, 1.0420 

II 0.3 0.0041,-1.0512 -0.0020, 1.0528 -0.0037, 1.0537 -0.0076, 1.0530 -0.0101, 1.0514 -0.0089, 1.0503 

I 0.1 -0.0486, -0.9975 0.0505, 0.9994 0.0486, 1.0008 0.0446, 1.0004 0.0421, 0.9985 0.0436, 0.9970 

II 0.1 0.0311, -1.0774 -0.0290, 1.0789 -0.0306, 1.0796 -0.0345, 1.0787 -0.0370, 1.0772 -0.0358, 1.0763 

I 0.06 -0.3913, -0.6543 0.3900, 0.6597 0.3816, 0.6681 0.3738, 0.6718 0.3736, 0.6676 0.3813, 0.6592 

II 0.06 0.0669, -1.1119 -0.0608, 1.1132 -0.0634, 1.1136 -0.0702, 1.1125 -0.0727, 1.1111 -0.0716, 1.1106 
II 0.04 0.1064, -1.1494 -0.1043, 1.1505 -0.1058, 1.1504 -0.1097, 1.1492 -0.1122, 1.1479 -0.1111, 1.1478 
II 0.02 0.1747, -1.2115 -0.1728, 1.2122 -0.1748, 1.2110 -0.1789, 1.2091 -0.1813, 1.2082 -0.1797, 1.2091 
II 0.005 0.2601, -1.2620 -0.2605, 1.2602 -0.2673, 1.2538 -0.2744, 1.2486 -0.2750, 1.2496 -0.2687, 1.2558 
II 0.001 0.3830, -1.1805 0.4036, 1.1593 -0.4828, 1.0840 -0.4252, 1.1420 

TABLE IV: Panorama of solutions as in Table but with sf 3 =0.01. Notice that the predicted masses, mi and m 2 have not 
been given since they are the same as in Tabic [I] 



type 


(m^min {S = 0) in 


cV (m 3 ) min (8 = 


0.5) (m 3 ) m4n (8 = 


1.0) (m 3 ) min (8 = 


1.5) (m 3 ) min (8 


= 2.0) (m 3 ) mi „ (S = 2.5) 


I 


0.0592716 


0.0592263 


0.0591312 


0.0590735 


0.0591074 


0.0592009 


II 


0.0006811 


0.0007755 


0.0009761 


0.0010993 


0.0010268 


0.0008287 



TABLE V: Minimum allowed value of the input mass, m 3 as a function of the input CP violation phase, 8 as in Table El but 



wherein, 

C = m\ (ci 2 c 13 ) 2 + m 2 (S12C13) 2 cos [2 (r 2 - n)] + m 3 (s 13 ) 2 cos [2 (r 3 - n + 6)] , 

D = m 2 (si 2 ci 3 ) 2 sin[2 (t 2 - ri)] + m 3 (si 3 ) 2 sin [2 (r 3 - ti + <5)] . (32) 

The dependence of \m ee \ on the input values of m 3 and 5, obtained by using the Ansatz of present interest, is displayed 
m Table. EHfor the same choices as in Table H] There is noticeable dependence on the input CP phase 5 for the larger 
values of ra 3 . 



type mi,m 2 ,m 3 in eV 


| m ee |(5 


= 0) m ce |(<5 


= 0.5) \m ee \(8 = 


1.0) |m ee |(<5 = 


1.5) |m ee |(<5 


= 2.0) |m ee |(<5 = 2.5) 


I 0.2955, 0.2956, 0.3 


0.164 


0.174 


0.183 


0.181 


0.170 


0.162 


II 0.3042, 0.3043, 0.3 


0.167 


0.177 


0.185 


0.183 


0.172 


0.164 


I 0.0856, 0.0860, 0.1 


0.051 


0.055 


0.057 


0.057 


0.054 


0.051 


II 0.1119, 0.1123, 0.1 


0.058 


0.062 


0.065 


0.064 


0.060 


0.057 


I 0.0305, 0.0316, 0.06 


0.026 


0.028 


0.029 


0.030 


0.029 


0.027 


II 0.0783, 0.0787, 0.06 


0.038 


0.040 


0.042 


0.041 


0.039 


0.037 


II 0.0643, 0.0648, 0.04 


0.029 


0.031 


0.032 


0.031 


0.029 


0.028 


II 0.0541, 0.0548, 0.02 


0.022 


0.023 


0.023 


0.023 


0.022 


0.021 


II 0.0506, 0.0512, 0.005 


0.019 


0.020 


0.020 


0.019 


0.019 


0.019 


II 0.0503, 0.0510, 0.001 


0.019 













TABLE VI: The neutrinoless double beta decay amplitude factor, \m ee \ in eV as a function of the input CP violation phase, 8. 
Here, the choice sf 3 = 0.04 has been made. 

For the reflected triangle solutions discussed above, the predictions of |m ee | are given below in Table EU A gam 
there is a noticeable dependence on S for the larger values of m 3 . However, the peak values occur at different values 
of S compared to Table VI. 

The effects of lowering sf 3 to 0.01 are finally illustrated below, for the non-reflected triangle case, in Table IVTTT1 
The main conclusion of this model for neutrinoless double beta decay, obtained by looking at all three tables above 

and noting the smooth dependence of \m ee \ on the inputs m 3 and S for each of the type I and type II solutions, is 

that \m ee \ should satisfy the restrictive approximate bounds: 



0.020eF < |m ee | < 0.185eV. 



(33) 
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type T7i3 in eV |m ee |(<5 = 0) \m ee \(S — 0.5) \m ee \(S — 1.0) \m ee \(S = 1.5) |m ee j(5 = 2.0) \m ee \(S = 2.5) 



I 


0.3 


0.164 


0.161 


0.167 


0.178 


0.183 


0.177 


II 


0.3 


0.167 


0.163 


0.169 


0.180 


0.186 


0.180 


I 


0.1 


0.051 


0.051 


0.053 


0.056 


0.058 


0.056 


II 


0.1 


0.058 


0.057 


0.059 


0.063 


0.065 


0.063 


I 


0.06 


0.026 


0.027 


0.029 


0.030 


0.030 


0.028 


II 


0.06 


0.038 


0.037 


0.038 


0.040 


0.042 


0.041 


II 


0.04 


0.029 


0.028 


0.029 


0.030 


0.032 


0.031 


II 


0.02 


0.022 


0.021 


0.022 


0.022 


0.023 


0.023 


II 


0.005 


0.019 


0.019 


0.019 


0.019 


0.019 


0.020 


II 


0.001 


0.019 













TABLE VII: The neutrinoless double beta decay amplitude factor \m ee \ in eV as a function of the input CP violation phase, 
5, using the reflected triangle. Notice that the predicted masses, 777,1 and m-2. have not been given since they are the same as in 
Table EH 



type 777,3 in eV |m ee |(5 = 0) |m ee |(<5 = 0.5) |m ee |(5 = 1.0) |m ee |(5 = 1.5) |m ee |(5 = 2.0) |m ee |(<5 = 2.5) 



I 


0.3 


0.177 


0.180 


0.181 


0.182 


0.179 


0.176 


II 


0.3 


0.179 


0.182 


0.184 


0.183 


0.181 


0.179 


I 


0.1 


0.056 


0.057 


0.057 


0.057 


0.056 


0.056 


II 


0.1 


0.063 


0.064 


0.064 


0.064 


0.063 


0.062 


I 


0.06 


0.029 


0.029 


0.030 


0.030 


0.030 


0.029 


II 


0.06 


0.041 


0.041 


0.042 


0.041 


0.041 


0.040 


II 


0.04 


0.031 


0.031 


0.031 


0.031 


0.031 


0.031 


II 


0.02 


0.023 


0.023 


0.023 


0.023 


0.023 


0.023 


II 


0.005 


0.020 


0.020 


0.020 


0.020 


0.020 


0.020 


II 


0.001 


0.020 


0.020 


0.020 






0.020 



TABLE VIII: The neutrinoless double beta decay amplitude factor |m ee | in eV as a function of the input CP violation phase, 
8, as in Table IvTl but with S13 = 0.01. Notice that the predicted masses, mi and 7712 have not been given since they are the 
same as in Table IVT1 



Here, the lower bound is intrinsic to the model but the upper bound reflects the experimental bound on the sum of 
neutrino masses quoted in Eq. Ij26(l and might be improved upon. We also note that, when both type I and type II 
solutions exist for a given value of 771,3, the type II solution gives somewhat larger |m ee |. The main dependence of 
\m ee \ is, of course, on the input parameter 7713. 

We plan to discuss elsewhere other physical applications including beta decay end point spectra and leptogenesis 
[23| using the approach of ref. [l5|. 



V. SYMMETRICAL PARAMETERIZATION 



In the parameterization of the leptonic mixing matrix given by Eq.©, which is similar to the one usually adopted, 
there appears to be an important distinction between the phase 5 and the phases Ti in the sense that only S survives 
when one considers the ordinary (overall lepton number conserving) neutrino oscillation experiments. However, this 
distinction may be preserved in a different way while using a more symmetrical parameterization. Thus consider, 



Ks = W23(#23, 023)^13(^13, 013)^12 (#12, 012), 

which contains the same mixing angles, as before but now has the three associated CP violation phases, 
ajij(9ij, 4>ij) were defined as in Eq. 10 • Writing out the whole matrix yields, 



K s = 



-Sl2C23e 



— *012 



C12C13 

Cl2Sl3S23e 



i(4>23— 013) 



Sl2Cl3e 



i<t>l2 



C12C23 - Sl2Sl3S23e 



»(012+023-013) 



si3e 



4<Pl3 



Si2S2 3 e^ (fe+fe) - c 12 s 13 c 23 e'^ ~c 12 s 23 e~^ - s^s^e^ 



Cl3S23e 

C13C23 



i<p23 



To relate this form to the previous one, we can use the identity 0, 

Ll>^ 1 (t)K S Ll>q(t) = ^23(6*23, ~ T 2 + 023 + T 3 )uJi 3 (9 13 , ~T\ + 0i 3 + T 3 )u>i 2 (012 , ~T\ + 012 + T 2 ), 



(34) 
The 

(35) 
(36) 
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where the diagonal matrix of phases, wq(t) was defined in Eq. (JJJ. Now choose the Tj's (two are independent) so that 
the transformed (12) and (23) phases vanish. Then, with the identifications, 



'12 


= n 


- T 2 


'23 


= T 2 


- 73 


'13 


= n 


- 7-3 



we notice that K$ is related to K of Eq. 



(37) 



K s =ujo(t)K. (38) 



Since K sits in the weak interaction Lagrangian, Eq. Q with the charged lepton field row vector on its left, all 
physical results will be unchanged if K is multiplied by a diagonal matrix of phases on its left. Thus K and Kg are 
equivalent and the relations between the CP violation phases of the two parameterizatons are given by Eq. (|37|l . By 
construction, S is the only CP violation phase which can appear in the description of ordinary neutrino oscillations. 
Solving the Eqs. I|37|) for 5 in terms of the three 0jj's gives, 



^13, 



(39) 



which shows that in the symmetrical parameterization, the "invariant phase" [24j combination I123 = 4>\2 + </>23 — </>i3 
is the object which measures CP violation for ordinary neutrino oscillations. It has the desired property of intrinsically 
spanning three generations, as needed for CP violation in ordinary neutrino oscillation or in the quark mixing analog. 
Furthermore, it can be seen [25| to have an interesting mathematical structure and to be useful for extension to the 
case where there are more than three generations of fermions 

The convenience of the symmetrical parameterization can already be seen in Eqs. Q32JI needed for obtaining \m ee \. 
These equations simplify when one observes that the combinations of phases occurring within them are simply 4>i2 
and 013. These are evidently the two phases which describe the coupling of the first lepton generation. Similarly, 
the estimates of CP violation needed for the treatment of leptogenesis made in section V of Jlj3] also simplify when 
expressed in terms of the 4>i j ■ 



VI. SUMMARY AND DISCUSSION 



Assuming that the squared neutrino mass differences and the three (CP conserving) lepton mixing angles are 
known, the mass of one neutrino and three CP violation phases remain to be determined. The complementary 
Ansatz, Tr{M u ) =0, provides two real conditions on the parameters of the lepton system. Here, we took and 5 as 
input parameters and then determined the other two masses and the other two CP violation phases according to this 
Ansatz. A geometric algorithm was presented based on a modification of an earlier treatment [l5T | in which the Dirac 
CP violation phase S was neglected, but the other two Majorana type CP violation phases were retained. That is a 
reasonable first approximation because the effect of delta is always suppressed by S13 which is known experimentally 
to be small. However, there is great interest in the determination of S so it should not be ignored. Additionally in ref. 
[THj l it was suggested that small CP violation scenarios might be close to the physical case. The present algorithm is 
exact and does not require the assumption that any parameters are small. 

The Ansatz yields a characteristic pattern for the neutrino masses and the CP violation phases. Because sf 3 is 
small, the main cause of change is the assumed value of the input parameter 7713. First one notes that the small 
experimental value of A in Eq. (|13fl always forces the neutrino 1 and neutrino 2 masses to be almost degenerate 
(See table HJ). For the largest allowed (from the cosmological bound Eq. 1)26(1 ) value of 7713, around 0.3 eV, there 
is an approximate three-fold degeneracy of all the neutrino masses. This is understandable since when the mass 
scale becomes large, both A and B can be considered negligible. Then the triangle of Fig. 1 becomes approximately 
equilateral. The internal angles of the triangle approximately measure the strength of the CP phases Ti — Tj and are 
clearly large in this situation. As m 3 decreases, a point around 0.06 eV is reached at which the type I solutions (7713 
largest) no longer exist. At this point the CP violation vanishes for the S =0 case and becomes small when S ^0. 
Also at this point the almost degenerate neutrino 1 and neutrino 2 masses decrease to about half the neutrino 3 mass 
in the type I case. In the present model neutrinos 1 and 2 never go below about half the mass of neutrino 3. The 
situation is a little different for the type II cases (7713 smallest). The type II solutions exist from a maximum of 7713 
about 0.3 eV to a minimum of about 0.001 eV. At the minimum the CP violation ceases for S = and has small 
effects when 5^0. Furthermore, at the minimum neutrinos 1 and 2 are about 50 times heavier than neutrino 3. 
Thus a possible hierarchy of neutrino masses can only exist in one way for the present model, with 777.3 considerably 
smaller than the other two. 
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Some further technical details of the model were displayed in Tables II -IV. These describe the S dependence of the 
limiting values of the input parameter m 3 just mentioned, the solution corrresponding to a reflected triangle and the 
effect of varying sf 3 . 

The model might be handy for getting an idea about the range of predictions for various leptonic phenomena, 
since it gives a plausible two parameter complete set of neutrino masses, mixing angles and CP violation phases. 
For example, in section V an application is made to the quantity \m ee \ which characterizes ncutrinoless double beta 
decay. The results for about the largest allowed sf 3 , the reflected solution and a smaller s^ 3 choice are shown in Tables 
VI, VII and VIII. It may be noted that the solutions vary rather smoothly with m 3 and 6 for a solution of given 
type. Thus, even though one might initially expect the result of allowing a two parameter family choice to be rather 
weak, it turns out that one gets fairly restrictive upper and lower bounds on \m ee \ as expressed in Eq. These 
approximate bounds may be considered also as a test of the present model. Of course, any direct determination of a 
neutrino mass, say from a beta decay endpoint experiment, will also provide a test of the model. 

In section V we took up a question which is independent of the present Ansatz. How should one parameterize the 
lepton mixing matrix? Of course, this is fundamentally a question of choice. However, we notice in the present work 
that the physical quantities we calculate, depend in the simplest way not on the conventional phases 5 and ti (where 
^2 T~i = 0) but on the quantities <pij given in Eq. i|37|) ■ Such a dependence arises naturally if one uses the alternative 
mixing matrix parameterization given by K$ in Eq. (|34|l . This may be understood physically in the following way 
The </>ij's by definition (see Eq.@ for example) span two generations. It is known [10( that CP violation begins at 
the two generation level for Majorana neutrinos. Thus it seems appropriate that the 4*^$, should appear. However 
if 8 = 0, the three (fiij's are not linearly independent according to Eq. (|39|) . This takes care of the two Majorana 
phases. When 5 ^ the three ifiij's are of course independent and the "invariant phase" combination /123 discussed 
some time ago 0,123 intrinsically spans three generations, as expected for "Dirac" type CP violation. 
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